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Abstract 

In  internally  manifolded  fuel-cell  stacks,  there  is  a  non-uniform  gas  flow  distribution  along  the  height  of  the  system  To  gain  an  insight  into 
this  distribution  an  analytical  model  has  been  developed.  In  the  model,  the  stack  is  viewed  as  a  network  of  hydraulic  resistances.  Some  of 
these  resistances  are  constant,  while  some  depend  upon  the  gas  velocity  and  can  be  determined  from  the  literature.  The  model  consists  of 
equations  for  the  network  with  counter-current  flow  in  the  manifold  channels.  Only  the  most  important  resistances  are  included,  i.e.,  the 
resistances  due  to  splitting  and  combining  the  flows  in  the  manifold  channels,  and  the  resistance  in  the  gas  channels  of  the  active  cell  area. 
The  ratio  between  the  average  flow  and  the  flow  in  the  upper  cell  can  be  solved  from  the  model.  In  this  manner,  a  very  useful  tool  for  separator- 
plate  design  is  obtained. 
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1.  Introduction 

Fuel  cells  produce  electricity  through  the  electrochemical 
conversion  of  gases.  Because  of  the  high  efficiency  and  the 
low  emission  level  of  pollutants,  fuel  cells  are  considered  to 
be  an  environment-friendly  way  of  producing  electricity. 

Depending  on  the  load,  one  single  fuel  cell  produces  0.5 
to  1  V.  To  yield  a  sufficiently  high  voltage,  the  cells  are 
stacked.  This  implies  that,  electrically,  they  are  connected  in 
series.  Stacking  of  the  cells  imposes  some  difficult  techno¬ 
logical  problems.  For  instance,  reactant  gases  have  to  be  able 
to  flow  on  either  side  of  each  cell.  This  requires  some  spacing 
but  that  the  electrical  contact  between  each  cell  has  to  remain. 
Also,  gases  may  not  leak  to  the  surroundings,  or  mix  with  the 
opposing  gas.  These  functions  are  combined  in  what  is  usu¬ 
ally  referred  to  as  the  separator  plate  or  interconnect. 

In  many  cases,  the  separator  plate  also  has  a  function  in 
the  manifolding  of  the  feed  gas  and  the  collection  of  the  waste 
gas.  This  is  the  case  of  internal  manifolding,  as  opposed  to 
external  manifolding  [1],  With  external  manifolding,  a 
header  is  placed  against  the  stack  of  cells  and  separator  plates. 
Both  these  technical  solutions  have  advantages  and  disad¬ 
vantages,  the  balance  of  which  still  has  to  be  determined. 

One  of  the  aspects  of  the  technology  of  stacking  of  fuel 
cells  is  how  well  the  gases  distribute  along  a  stack  of  fuel 
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cells.  A  non-uniform  distribution  will  <  ause  differences  in 
the  performance  of  each  cell,  which  is  allowed  within  certain 
limits  [  1  ] .  Important  parameters  in  such  a  prediction  are  the 
actual  gas  flows  and  the  resistance  that  the  gas  meets  when 
flowing  through  the  system. 

The  model  presented  in  this  paper  provides  an  analytical 
approximation  and  is  applicable  for  internally  manifolded 
fuel-cell  stacks.  The  feed  gases  and  waste  gases  in  the  man¬ 
ifold  channels  are  assumed  to  flow  in  opposite  directions, 
which  technically  is  the  most  viable  option. 

The  analytical  approach  has  the  disadvantage  that  simpli¬ 
fications  need  to  be  made  to  get  a  result.  The  advantage  is 
that  under  the  assumptions  (or  simplifications),  the  result  is 
verifiable,  as  opposed  to  a  numerical  analysis.  Furthermore, 
the  analytical  approach  clearly  identifies  the  importance  of 
each  parameter  in  relation  to  the  gas  distribution. 


2.  Theory 

The  gas  flow  distribution  in  a  fuel  cell  stack  is  determined 
by  the  hydraulic  resistances  in  the  separator  plate  and  mani¬ 
fold  channels.  The  most  important  resistance  in  the  separator 
plate  is  caused  by  the  channels  adjacent  to  the  active  cell  area. 
Other  important  resistances  are: 
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(i)  the  resistance  caused  by  the  flow  through  the  manifold 
channels;  the  manifold  channels  can  be  regarded  as  pipes 
with  a  rough  surface; 

(ii)  the  resistance  that  results  from  the  pressure  drop 
caused  by  splitting  of  the  flow  in  the  inlet  manifold  channels, 
and 

(iii)  the  resistance  that  results  from  the  pressure  drop 
caused  by  combining  the  flow  in  the  outlet  manifold  channels. 

The  resistance  in  the  channels  adjacent  to  the  cell  area  can 
be  inferred  from  laminar  flow  formulae  [2,5].  Flow  resis¬ 
tance  in  pipes  with  a  rough  surface  is  described  in  Refs.  [2- 
4  ] ,  while  resistance  data  for  splitting  and  combining  are  given 
in  Refs.  [2,4] .  The  resistance  depends  upon  the  fraction  that 
is  split  off  or  added  to  the  mainstream  in  the  manifold 
channels. 

Flow  results  from  a  pressure  drop.  For  a  simple  one-dimen¬ 
sional  duct  flow,  the  relation  between  flow  and  resistance  can 
be  given  by 


2.2.  The  resistance  in  the  manifold  channels 

The  pressure  drop  between  two  ad  jacent  separator  plates 
in  the  inlet  and  outlet  manifold  channels,  A and  A pmo 
respectively,  can  be  written  as 

A  =  (4) 

Apmo  —  K2<f>l  (5) 

where  and  <f>0  are  the  flows  in  the  manifold  channels  of  the 

inlet  and  outlet,  respectively. 

Kx  represents  the  sum  of  the  resistance  of  the  rough  pipe 
flow  and  the  resistance  due  to  splitting,  while  K2  represents 
the  sum  of  the  resistance  of  the  rough  pipe  flow  and  the 
resistance  due  to  combining,  and  are  described  by  Eqs.  (6) 
and  (7): 


where  Ap  is  the  total  pressure  drop  \  K  the  hydraulic  resis¬ 
tance  coefficient,  p  the  density,  v  the  flow  velocity  averaged 
over  the  cross-sectional  area,  /  the  flow  length,  and  Dh  is  the 
hydraulic  diameter  of  the  channel. 

2. 7.  The  resistance  in  the  cell  channels 


The  flow  in  these  channels  is  assumed  to  be  laminar.  This 
is  valid  for  the  anode  and  cathode  gasflow  of  an  adiabatic 
solid  oxide  fuel  cell  (SOFC)  with  a  flow  length  of  up  to 
1  m,  and  can  be  checked  by  calculating  the  actual  Reynolds 
number  for  a  specific  case. 

For  laminar  flow,  K  is  inversely  proportional  to  v  and, 
therefore,  the  pressure  drop  in  Eq.  (1)  is  linear  with  v  [2], 
and  hence  it  can  be  shown  that 


Apc  =  Kc<f> 


with 

ic 

c  kAQ  Dl 


(2) 

(3) 


where  17  is  the  gas  viscosity,  k  the  number  of  channels  in 
parallel  that  can  be  allocated  to  one  manifold  channel,  lc  the 
corrugation  length,  Ac  the  cross-sectional  area  of  one  corru¬ 
gation,  Dhc  the  hydraulic  diameter  of  one  corrugation,  and  <f> 
the  flow  that  passes  the  cell  area  considered  (which  is  equal 
to  kvAc) ;  /  is  a  friction  factor  that  depends  on  the  channel 
geometry.  For  a  channel  with  a  circular  cross  section  f—  32, 
for  a  channel  with  a  triangular  cross  section/=  26.67,  and  for 
a  channel  with  a  square  cross  section/=  28.4  [5] . 

1  The  sum  of  static,  dynamic  and  gravitational  pressures. 


where  K ^  and  Kmo  are  the  resistance  coefficients  for  splitting 
and  combination  of  the  flow,  respectively,  and  Kp  is  the  resis¬ 
tance  coefficient  for  the  rough  pipe  flow  of  the  manifold 
channel.  These  coefficients  are  dimensionless  and  are  tabu¬ 
lated  in  handbooks  [2,4]  or  have  to  be  inferred  from  meas¬ 
urements.  Dhm  is  the  hydraulic  diameter  of  the  manifold 
channel,  /m  is  the  distance  between  two  separator  plates, 
measured  between  the  centres  of  the  plates,  and  Am  is  the  area 
of  a  manifold  channel.  It  is  assumed  that  the  inlet  and  outlet 
diameters  are  equal. 

23.  Model 

The  stack  can  be  represented  as  a  network  of  hydraulic 
resistances,  in  parallel  and  in  series,  as  depicted  in  Fig.  1 .  In 
this  figure,  a  subscript  is  added  to  the  symbols  used  for  the 
flow,  to  indicate  the  cell  number.  The  convention  for  the 
symbols  for  pressure  is  analogous.  It  is  assumed  that  all  sep¬ 
arator  plates  are  of  equal  design  and  geometry. 

Using  Eqs.  (2),  (4)  and  (5),  the  pressure  drop  across  the 
cell,  with  index  j,  can  be  written  as 

Pij-Py-Kt{  (4>j+  ■■■  +  <£„)2+  (<£,+  1  + ...  +  <t>n)2 
+  ...  +  (<t>n)2}+Kc<f>n  +  KM<j>n)2 

+  (0/.  +  0n-t)2 +•••  +  (<£..  +  + 

=  (AT,  +*2>E<i>)2  +  *A  (8) 

m—j  i  —  m 

The  terms  between  the  round  brackets  in  the  first  line  of 
Eq.  (8)  equal  the  local  flow  in  the  manifold  inlet  channel. 
The  second  line  is  similar  but  relates  to  the  outlet  channel. 

In  Eq.  ( 8 )  it  is  already  assumed  that  and  K2  are  constant. 

Actually  the  parameters  depend  on  the  fraction  that  is  split 
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Fig.  1 .  Flow  through  a  stack  represented  in  a  network.  Kx  and  K2  are  the  inlet 
and  outlet  manifold  resistances;  Kc  is  the  resistance  of  the  corrugations; 
is  the  flow  through  the  channels  of  cell  j\  pXj  and  p2J  are  the  respective 
pressures  at  the  in*  and  outlet  of  cell  j. 


Fig.  2.  Flow  distribution  for  stacks  with  countercurrent  manifold  flow. 

from,  or  added  to,  the  mainstream.  This  fraction  is  small, 
however,  and  therefore  approximately  constant  for  all  cells, 
except  for  the  upper  cells.  Therefore,  constant  resistance  coef¬ 
ficients  are  assumed.  By  doing  this,  most  of  the  inaccuracy  is 
introduced  in  the  top  of  the  stack. 

To  approximate  the  sum  of  the  terms  in  Eq.  (8),  it  is 
assumed  that  the  distribution  is  rather  good.  This  is  done  by 
substituting  <f}t  =  </>,  which  is  the  average  flow  in  cells  i  ton. 
Eq.  (8)  will  then  become 

Pm -Pv = Kc <t>n  -  ( *,  +  jra)£( 2/  - . . . ) 2  (9) 

O 

From  the  network  in  Fig.  2,  it  follows  that 

2  For  brevity,  read  for  ther  last  term:  [  2/  —  3/2(2n  4-  3 )  +  j ( 6n2  +  18n  + 
13)  —  (/i+  1  )(2n2  +  7rt +  6)]. 


Pm~Pv=kc 

Thus,  combining  Eqs.  (9)  and  ( 10)  gives 


Kc  6 


(10) 


(11) 


Assume  that  the  average  flow  for  the  entire  stack  (i.e.,  for 
i  =  1  to  n)  is  found  in  cell  J,  then  it  can  be  written  as 


Substituting  Eq.  (11)  results  in 


(12) 


(13) 


By  comparing  Eq.  (11),  for  J=j  —  1,  and  Eq.  (13),  the  cell 
that  gets  the  average  flow  is  found.  Table  1  shows  the  number 
of  plates  and  the  relative  location  of  the  cell  with  the  average 
flow.  For  stack  heights  of  interest  (n  >  1 00),  this  is  the  case 
at  approximately  38%  of  the  stack  heigh  i. 

It  is  also  now  known  that  cell  J  gets  the  average  flow.  This 
information  is  required  in  order  to  determine  the  following 
flow  difference  ratio 


r.fcA.(2(/41),-.)  (14) 

c pi  —  efrn  n(2n2  —  3n+l) 

This  ratio  gives  the  quotient  of  the  flow  difference  between 
the  average  and  the  top  cell,  and  the  bottom  and  the  top  cell. 

Substituting  7/n,  as  calculated  in  Table  1,  gives  the  ratio 
as  a  function  of  n.  In  Table  2,  the  values  of  r  are  listed  for 
several  values  of  n.  The  latter  data  show  that,  for  practical 
purposes,  it  may  be  assumed  that  r  =  0.25.  This  indicates  that 
25%  of  the  cells  receive  less  than  the  average  flow,  while  the 
remainder  receive  more. 


Table  1 

Relative  position  of  the  cell  with  the  average  flow  lor  various  stack  heights, 
n 


n 

Jfn 

20 

0.395 

50 

0.380 

100 

0.375 

200 

0.372 

1000 

0.370 

Table  2 

Ratio,  r,  for  various  stack  heights,  n.  The  ratio,  t  indicates  the  fraction  of 
the  cells  that  receive  less  than  the  average  flow  in  i  he  stack 


n 

r 

20 

0.289 

50 

0.257 

100 

0.247 

200 

0.243 

1000 

0.239 
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Eq.  (13)  can  be  rev.  ;itten  as  follows  (realise  that  </)j  =  <f>) 


~Ki± 
_  K 


+  K2  cf> 


12 


n(n2—  1) 


(15) 


The  parameters  on  the  right-hand  side  of  Eq.  (15)  follow 
from  the  geometry  of  the  separator  plate  and  the  number  of 
cells,  and  the  amount  of  gas  supplied.  The  quotient  <£n/<£ 
gives  the  ratio  between  the  flow  through  the  top  cell  and  the 
average  flow.  For  instance,  if  a  flow  is  set  for  an  average 
utilization  of  80%,  and  may  not  exceed  90%  utilization  in  the 
top  cell,  then  4>n/ 4>  becomes  80/90  =  0.89.  If  a  lower  value 
is  found  on  the  right-hand  side  of  Eq.  (15),  then  the  maxi¬ 
mum  utilization  is  more  than  90%.  This  means  that  the  num¬ 
ber  of  cells  must  be  decreased,  or  the  geometry  of  the 
separator  plate  must  be  adapted. 


3.  Case  study 

To  show  the  relevance  of  Eq.  ( 15)  a  case  study  has  been 
undertaken.  Unfortunately,  not  many  details  of  the  design  of 
separator  plates  are  available  in  the  open  literature.  Hence  it 
was  decided  to  use  what  is  available,  and  make  assumptions 
about  the  details  not  available.  A  good  basis  to  work  from 
was  found  in  Refs.  [6,7],  where  the  technology  of  Tokyo 
Gas  is  shown. 

Assume  that  an  average  utilization  of  80%  is  targeted  and 
that  it  is  not  wanted  to  exceed  90%  utilization  in  the  top  cell. 
As  indicated  above,  the  right-hand  side  of  Eq.  ( 15)  has  to  be 
greater  than  0.89,  which  restricts  the  number  of  cells,  n. 
Table  3  shows  the  data  used  for  the  calculations. 

The  gas  flow  (Table  3)  was  calculated  on  the  basis  of  an 
80%  average  utilization  (the  gas  was  pure  hydrogen),  the 
cell  surface  area  was  assumed  to  be  100  cm2,  and  a  current 
density  of  1  A/ cm2  was  obtained  from  the  cell,  i.e.,  the  current 
density  at  maximum  power  output  [7].  It  should  be  noted 


Table  3 

Data  used  for  calculation  of  critical  stack  height  for  the  anode  gases 


Symbol 

Name 

Value 

Unit 

k 

channels  /  manifold 

50 

_ 

lc 

channel  length 

0.1 

m 

Ac 

channel  cross  section 

1X10~6 

m2 

AlC 

channel  hydraulic  diameter 

1X10'3 

m 

/ 

friction  factor 

28.4 

— 

L 

separator  plate  pitch 

0.005 

m 

Aim 

manifold  hydraulic  diameter 

0.02 

m 

manifold  cross  section 

3.5  X10-4 

m2 

Xml 

resistsance  coefficient  for  splitting 

0.16 

— 

Ana 

resistance  coefficient  for 
combining 

0.16 

— 

xp 

resistance  coefficient  for  rough 
pipe  flow 

0.14 

<t> 

flow  per  cell  per  manifold,  at  T 
and  p 

6.77X10“ 5 

m3/s 

P 

average  gas  density,  at  T  and  p 

0.074 

kg/m3 

V 

viscosity,  at  T 

4.81  X10“5 

Pa  s 

that  the  flow  has  to  be  used  at  the  actual  temperature  and 
pressure.  The  current  density  appears  to  be  very  high,  but  in 
Ref.  [7]  it  is  indicated  that  this  is  the  direction  where  the 
technology  is  heading. 

To  calculate  the  pressure  drop  ac  ross  a  cell,  it  must  be 
realised  that  the  gas  viscosity  changes  due  to  the  electrochem¬ 
istry.  At  the  cell  inlet  is  hydrogen,  w  hile  at  the  cell  outlet  is 
a  mixture  of  hydrogen  and  water  vapour,  the  average  com¬ 
position  of  which  is  20%  hydrogen  and  80%  water  vapour. 
Also  depending  on  its  location  in  the  stack,  each  cell  has  a 
different  outlet  composition  and,  consequently,  a  different 
viscosity  and  a  different  value  of  Kc  For  the  calculation,  the 
average  composition  and  the  average  mix  viscosity  of  inlet 
and  outlet  gases  were  used. 

Taking  the  values  in  Table  3,  and  the  above  assumptions, 
the  right-hand  side  of  Eq.  (15)  reaches  the  critical  value  at  a 
stack  height  of  76  cells.  At  this  stack  height,  the  top  cell  has 
a  utilization  of  90%.  The  cell  that  gets  the  average  flow 
(which  is  at  approximately  38%  of  the  stack  height)  is  cell 
no.  29,  counted  from  the  fuel  gas  supply  side. 

With  the  ratio  r,  it  is  found  that  the  flow  in  the  bottom  cell 
is  1.73  times  the  average  flow.  This  implies  that  the  flow 
through  the  bottom  cell  is  almost  twice  the  flow  through  the 
top  cell  for  this  number  of  cells. 

A  similar  calculation  was  made  for  the  cathode.  Assuming 
adiabatic  operation,  an  average  temperature  increase  of 
150  °C,  and  the  same  critical  factor  as  for  the  anode,  it  is 
found  that  the  critical  height  is  reached  at  18  cells.  The  high 
flow  required  for  adiabatic  operation  (almost  20  times  stoi¬ 
chiometry)  is  the  main  cause  for  the  lower  cell  number. 

The  results  of  the  calculations  depend  strongly  on  the  val¬ 
ues  of  Kpy  Km j  and  Kmo.  The  values  used  here  are  tabulated  in 
handbooks  [2-5]  and  are  related  to  a  geometry  that  differs 
from  the  actual  stack  geometry.  By  actually  measuring  the 
pressure  drops  in  a  stack,  however,  more  accurate  values  will 
be  determined.  These  measurements  can  also  be  carried  out 
at  room  temperature,  ambient  pressure,  and  using  air  in  a 
stack  with  dummy  cells. 


4.  Conclusions 

1.  The  flow  distribution  for  fuel  cell  stacks  can  be  given  in 
analytical  form,  thereby  a  very  useful  tool  is  obtained  for 
separator-plate  design. 

2.  The  cell  that  is  at  38%  of  the  stack  height  receives  a  flow 
that  is  equal  to  the  average  flow. 

3.  About  25%  of  the  cells  receive  less  than  the  average 
flow,  the  rest  receive  more. 


5.  List  of  symbols 

Ac  cross-sectional  area  of  one  gas  channel  adjacent  to 
a  cell,  m2 

Am  cross-sectional  area  of  one  manifold  channel,  m2 
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Dhc  hydraulic  diameter  of  one  gas  channel  adjacent  to 
a  cell,  m 

Aim  hydraulic  diameter  of  one  manifold  channel,  m 
/  friction  factor,  — 

k  number  of  channels  allocated  to  a  manifold,  - 

Kx  resistance  factor  relating  the  effect  of  rough  pipe 
flow  and  splitting  of  a  flow  to  a  pressure  drop, 
Pa/(m3/s)2 

K2  resistance  factor  relating  the  effect  of  rough  pipe 
flow  and  combining  of  flows  to  a  pressure  drop, 
Pa/(m3/s)2 

Kc  resistance  factor  relating  laminar  flow  to  pressure 
drop,  Pa/(m3/s) 

resistance  coefficient  resulting  from  splitting  a 
flow  in  two  flows,  — 

Kmo  resistance  coefficient  resulting  from  combining 
two  flows  in  one,  — 

Kp  resistance  coefficient  for  rough  pipe  flow,  - 
lc  length  of  gas  channel  adjacent  to  a  cell,  m 

separator  plate  pitch,  i.e.,  distance  between  two 
adjacent  separator  plates,  measured  between  the 
plate  centres,  m 

n  number  of  cells  in  a  stack,  - 

p  pressure,  Pa 

v  gas  velocity,  m/s 

Greek  letters 

<t>  flow  per  cell  per  manifold  at  T  and  p,  m3/s 
7)  viscosity  at  T,  Pa  s 

p  average  gas  density  at  T  and  p,  kg/m3 


Subscripts 

1  relating  to  the  inlet  manifold  c  hannel 

2  relating  to  the  outlet  manifold  channel 

j  cell  number  in  a  stack  counted  from  the  gas  supply 

side 

J  cell  number  in  a  stack  counted  from  the  gas  supply 
side  that  gets  the  average  flow 
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